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Abstract 

<N 

Almost all organisms show some kind of time periodicity in their behavior. Especially in mammals 
the neurons of the suprachiasmatic nucleus form a biological clock regulating the activity-inactivity 
cycle of the animal. This clock is stimulated by the natural 24-hour light-dark cycle. In our model of 
this system we consider each neuron as a so called phase oscillator, coupled to other neurons for which 
the light-dark cycle is a Zeitgeber. To simplify the model we first take an externally stimulated single 
t-H phase oscillator. The first part of the phase interval is called the active state and the remaining part 

is the inactive state. Without external stimulus the oscillator oscillates with its intrinsic period. An 
external stimulus, be it from activity of neighboring cells or the periodic daylight cycle, acts twofold, 
it may delay the change form active to inactive and it may advance the return to the active state. 
The amount of delay and advance depends on the strength of the stimulus. We use a circle map as 
a mathematical model for this system. This map depends on several parameters, among which the 
intrinsic period and phase delay and advance. In parameter space we find Arnol'd tongues where the 
system is in resonance with the Zeitgeber. Thus already in this simplified system we find entrainment 
i i and synchronization. Also some other phenomena from biological experiments and observations can 

be related to the dynamical behavior of the circle map. 

> 
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1 Introduction 

O Setting of the problem 

Rhythmic behavior is present in almost all organisms. Their rhythms can be autonomous, but are more 
often externally stimulated. One such stimulus is the 24-hour natural light-dark cycle which governs 
the activity-inactivity cycle of many animals and plants. The latter is the most common Zeitgeber or 
periodic stimulus, although an alternating high-low temperature cycle is another example of a Zeitgeber. 
Shaped by several millions of years of evolution under this 24-hour light-dark cycle |22j , many organisms 
still exhibit these behavioral rhythms (often with slightly different period) even in conditions without 
information on the alternation of light and darkness. How the almost-24-hour intrinsic period of the 
internal rhythm synchronizes with or entrains to an external Zeitgeber is one of the major questions in 
circadian biology, see [H [10] . 

In mammals, the circadian clock resides in the suprachiasmatic nucleus, a neuronal hypothalamic tissue 
residing just above the optic chiasm. It consists of about 10 000 interconnected neurons or pacer cells 
[T71 |2"T1 l2"rl , 29 . Experimental evidence seems to support the model of the suprachiasmatic nucleus as 
a collection of so called phase oscillators, which was put forward in 1980 to explain circadian rhythms 
[14j . Each oscillator has its own intrinsic period of activity-inactivity [36] and interactions among them 
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are believed to synchronize their oscillations [3U1 UHl 132 ■ This model is used to explain the observation 
that the activity-inactivity cycle of an organism closely follows the period of a Zeitgeber like the 24-hour 
light-dark cycle. 

However, before trying to analyze a model for the collection of pacer cells, we consider a model for a 
single cell. A collection of interacting pacer cells may mathematically be modelled by a (large) number of 
coupled almost identical oscillators. Such a model is rather involved, therefore we first make the following 
simplifying assumption. In our model each oscillator experiences the external (averaged) forcing by the 
collection of oscillators but does not influence the dynamics of the others. Thus the coupling we consider 
is asymmetric. This approach is also widely used in astro-dynamics . The basis for our model of a single 
pacer cell is given in [6l [14] where the state of a pacer cell is determined by the phase in its activity- 
inactivity cycle. The stimulus of an external Zeitgeber may advance or delay the phase depending on the 
phase itself but also on parameters characterizing the pacer cell. The parameters essentially consist of 
the intrinsic period of the pacer cell, the intrinsic length of the activity interval and the strength of the 
interaction with an external stimulus. The latter can be the relative number of active pacer cells in the 
environment as well as a stimulus originating from an external light, temperature, etc level. In our model 
the Zeitgeber will be a periodic, quasi-periodic or even more general function of time. But in this article 
we restrict to a periodic Zeitgeber where the period is an average or prevailing period in the natural light 
and dark cycle or in the collective behavior of pacer cells in the environment. 

Models based on phase oscillators date back to at least 1967, see [10], and have been studied by several 
others [TBI [28] . For an overview see [T5] . 

Main questions 

Thus we study a model of a single pacer cell, stimulated by its environment but not contributing to the 
collective behavior. For such a situation the main questions we wish to address are 

1. Can a single pacer cell synchronize with or entrain to a periodic Zeitgeber? 

2. If so, how does this depend on properties of the pacer cell? 

The mathematical model we use to describe a pacer cell stimulated by a periodic Zeitgeber, is a dynamical 
system, more specifically a map on the circle. This map depends on the parameters characterizing the 
pacer cell. Typical dynamics of a circle map most relevant in view of the questions above is dynamics 
of fixed points and dynamics of periodic points. Fixed points correspond to entrainment of the pacer 
cell which means that the sequence of onset times of the activity interval has the same period as the 
Zeitgeber. Periodic points correspond to synchronization meaning that there are an integer number of p 
different onset times of the activity interval during another integer number of q periods of the Zeitgeber. 
Such points are called p : q periodic points. In this vocabulary fixed points are 1 : 1 periodic points, in 
other words entrainment is a special kind of synchronization. Then the questions for the single pacer cell 
are translated into the following questions for the mathematical model. 

1. Do stable fixed or periodic points exist for the map on the circle? 

2. If so, for which domain in parameter space? 

Summary of results 

The analysis of the model for a single pacer cell shows that both entrainment and synchronization are 
possible. For weak interaction with the Zeitgeber the parameter space is divided in regions with unique 
and distinct dynamics. When we fix all parameters except the intrinsic period of the cell and the strength 
of the interaction we get the so called Arnol'd tongues, see figure [6| These are labelled by p and q such 
that for parameter values in the p : q tongue, the dynamics of the map is p : g-periodic corresponding to 
synchronization of the pacer cell. This shows that for a fixed strength of interaction with the Zeitgeber 
there are ranges of synchronization. Which means that for only a finite interval of values of the intrinsic 
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period synchronization is possible. A phenomenon which is confirmed by biological experiments and 
observations, see [21 IH IS IISl E21 HH [321 HH 133 ■ As a remark we note that varying the 

intrinsic period at a fixed period of the Zeitgeber is equivalent to varying the period of the Zeitgeber 
while the intrinsic period is fixed. The largest range is the range of entrainment since the 1 : 1 or main 
tongue is the largest. The location of the range of entrainment depends on the parameters governing 
phase advance and delay. In one extreme case only pacer cells with an intrinsic period smaller than the 
period of the Zeitgeber can be entrained, in the other extreme case only pacer cells with a larger intrinsic 
period can be entrained. Furthermore the onset time of the activity interval is an increasing function of 
the intrinsic period see section |4j which has been observed in many organisms, see [51 \2[ IST)) 155]. 

2 Biological model 

In our biological model of pacer cells, the state of each pacer cell is determined by a single variable: a 
phase 9 <E [0, r]. This phase is the resultant of many biochemical processes. In an isolated pacer cell the 
phase increases in time with speed one until it reaches a value r, then it jumps to zero and starts to 
increase again. If the phase is between zero and a value a < r the cell is active and between a and r, the 
cell is inactive. Thus an isolated pacer cell shows a periodic activity-inactivity cycle with period r. We 
call a the length of the intrinsic activity interval and r is called intrinsic period of the pacer cell. Both 
the values of a and r are properties of the individual cell. 

In a collection of pacer cells, considered as an organ without external stimulus, the interaction among the 
cells is modelled as follows. Let t n be the time at which the phase is zero. Then the cell is active during 
the time interval t n to t n + a. At the end of this interval an isolated cell would become inactive but now 
the cell remains active until a later time namely t n + a + sZ, where Z is the fraction of active cells and 
e is a small cell dependent parameter. We incorporate this phase delay into the model by once delaying 
the phase by an amount of eZ at time t n + a, so instantly 9 becomes 9 — eZ . Similarly the inactivity 
period of the cell is shortened by activity of other cells. If at time t n+ i the phase of the cell is below r 
by an amount of rjZ so that 9 + r/Z = t the cell becomes active again. Thus there is a phase advance of 
9 by an amount of r/Z at time t n+ \, that is instantly 9 becomes 9 + r/Z. Since the latter is equal to r, by 
previous assumptions the phase jumps to zero. 

Before isolating a single pacer cell we note that in a collection of n pacer cells the phase of cell i 
is 9i. The quantity Z for cell i depends on all phases 9k except 9i and possibly also on time t, so 
Z = Zi(6\, . . . , 9i, . . . , 9 n , t), where the hat means 9i excluded. However for reasons to be explained 
below, we will later on consider Z as a function of time only. 

In a model with a large number of interacting pacer cells we assume that the interaction is uni-directional. 
That is the fraction of active cells stimulates a single cell, but the influence of a single cell on the collection 
is negligible. Therefore the quantity Z can be considered as only time dependent Z — Z(t). Now our 
final model consists of a single pacer cell and we study its response to the external stimulus Z. The pacer 
cell is characterized by four parameters e, rj, a and r. We do not make a distinction whether Z is related 
to the fraction of active cells as a function of time or a daily light and dark signal. In our model Z is the 
Zeitgeber for the cell and we will discuss several choices for it. 

In principle the Zeitgeber Z can be any function of time. Here we will take it periodic, where the period 
is for example an average of the observed periods of the collection of pacer cells or an average period 
of an external stimulus like the daily light and dark cycle. But we may also take a much longer period 
to model seasonal effects. Since pacer cells may have different values of s, rj, a and r we are especially 
interested in the domain in parameter space where synchronization occurs. 



3 Mathematical model 

The system described in section [2] has a state which evolves in time. This calls for a dynamical systems 
approach, as a general reference see [9]. In order to do so we have to identify a state space and an 
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evolution law. This means that we should find one or more quantities describing the state of the system 
and then find an evolution law that uniquely defines a future state once an initial state is given. Here 
we first restrict to deterministic dynamical systems that is we do not include noise or some other kind 
of random fluctuations. It will turn out that we can define a discrete dynamical system on the circle. 
Using this model we are able to answer some of the questions in section [l] for a single pacer cell with a 
Zeitgeber. 

An isolated pacer cell. The model we build is based on the phase 9 of the cell. First we consider the 
cell as a 'free running' oscillator, depending on the two parameters a and r. Later on we also include a 
periodic Zeitgeber with two additional parameters e and rj governing the strength of the forcing. 

Let 9 : K — > [0, r) be the phase of the cell given by 




t - t n , for t G [t n ,t n+ x) 
0, forf = i n+1 , 



where r is the maximal value of 9. When 9 reaches the value a, with < a < r, the cell undergoes 
the transition from active to inactive. Suppose at t = t n the phase 9 is zero. The phase increases with 
speed 1 until it reaches the maximal value r, then 9 instantly becomes zero again. This happens at 
time t = t n+ i = t n + t. Thus the period of the 'free running' oscillator is r. The transition times are 
t n = to + nr and t n + a where the cell state changes from inactive to active and from active to inactive 
respectively, see figure [T] 




t n + Oi t n +T t n + a t n +T 



Figure 1: The phase 9 as a function oft without (left) and with (right) Zeitgeber. Suppose the phase 9 
is at time t n , then the state of the cell is active for < 9 < a and inactive for a < 9 < r . Recall 
that a + p = t. Without Zeitgeber (left figure) this corresponds to the time intervals [t n ,t n + a) and 
[t n + a, t n + r). In presence of a Zeitgeber the cell is active in the time interval [t n , t n + a + sZ(t n + a)) 
and inactive in the time interval [t n + a + eZ(t n + a), t n+ i), now t n+1 is defined implicitly, see text. 



We will use the transition times t n in phase space K to define a dynamical system whose evolution takes 
t n into t n+ i. Without Zeitgeber we have a sequence t n = to + nr as described above depending on the 
parameters a (although trivially) and r. Next we include a Zeitgeber. 

A single pacer cell with a periodic Zeitgeber. To model a non-isolated pacer cell in a collection of 
other pacer cells with or without an external stimulus, we consider a single pacer cell with a Zeitgeber Z 
which is a function of time only. First we define the Zeitgeber Z. Though non-essential it is convenient 
to scale time so that the period of the Zeitgeber becomes one. 

Definition 1. The positive function Z : M. —¥ [0, 1] satisfies the following 



i) Z is differentiate, 
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ii) Z is periodic with period 1. 

In the presence of a Zeitgeber the phase 6 again increases with speed 1, starting at 9 = at t = t n , but 
when 9 reaches the value a it instantly drops back by an amount of eZ(t n + a). Then it again increases 
with speed 1 until it reaches a value at t = t n +i such that 9(t n +i) + T)Z(t n +i) — r. Note that i n +i is 
implicitly defined. Thus the phase 9 of the cell is given by 

!t - t n , for t € [t n ,t n + a) 

t-t n -eZ(t n + a) for t e [t n + a,t n+1 ) (1) 
0, for * = t n+ i. 

Apart from the parameters a and r we now also have e and rj. The latter two give the 'strength' of the 
Zeitgeber. In order that the model be consistent we impose the following conditions on the parameters 

< a < t, £ > 0, 77 > 0, a — e > 0, a — £ + 77 < r, 

so that 9 remains between and r. 

Dynamical system. The state of the dynamical system we define is the transition time t n rather than 
the phase 9. Indeed solving equation 

9(t) + r\Z{t) = t or equivalently t - t n - eZ(t n + a) + rjZ(t) = t (2) 

for t, under conditions to be specified later, yields a unique solution t = t n+ \ once t n is given. Thus we 
may write t n+ i — F^(t n ) for a map : R — > R depending on parameters /1 = (s, r), a, r). However, it 
turns out that has the property F^(t+1) = F^(t) + 1 so that F^ is the lift of a circle map f^-.S^-^S 1 . 
This means that we now have a dynamical system with phase space the circle S 1 and evolution law f^. 
Conceptually it is easier to work with the circle map f ^ but for actual computations we usually prefer 
the lift Fp. Let us summarize the result in the following proposition, for a proof see the appendix. 

Proposition 2 (Circle map and lift). Let Z be as in definition^ and define the function U e : R — > R as 
U e (t) = t + eZ(t). Then the map F^ : R ->• R with 

F^t) = U-\U e {t + a)-a + T) (3) 

defines a parameter dependent differentiate dynamical system, provided that r\ is small enough. Param- 
eters are /1 = (e, 77, a, t). Furthermore, F^ is the lift of a circle map / M : S 1 — > S 1 of degree one, given by 



U(t) = U-\U E (t + a)- ol + t) modi. (4) 

A circle map is a one-dimensional map just like an interval map. The dynamics of the two have much in 
common which is most prominent when the circle map is studied by a lift. However, because the circle is 
different from the interval, there are also differences in dynamical behavior. For example non-degenerate 
fixed points of a circle map come in pairs. See figure [2| 

We now make a further distinction between two cases, namely whether /„ is invertible (a diffeomorphism) 
or not (an endomorphism) . The difference between these cases is not only in dynamical behavior but 
the second case has far richer bifurcation scenarios. Essentially it boils down to both e and 77 being 
'sufficiently small' or one of them not meeting this criterion. However, there is a priori no reason to 
assume that either of them is small, (nog meer over te zeggen?) 

a) F^ is the lift of a circle diffeomorphism. In this case F^ is differentiable and F^ 1 exists and 
is also differentiable. Both and U e have to be invertible, which in turn means that e and 77 must 
be small enough. 

b) F^ is the lift of a circle endomorphism. In this case F M is again differentiable but F^ 1 does 
not necessarily exist. Now only U n has to be invertible, which means that only rj must be small 
enough. 
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o t* s f* i 

Figure 2: Phase portrait of circle map and graph of lift F^. f ^ has two fixed points indicated by bullets. 
One is stable, the other is unstable according to the arrows. The lift F^ is a map on the interval [0, 1] 
with F^Y) = F^(0) + Y, drawn is Fp — Y. t* (stable) and £* (unstable) satisfy F^(t) = t + Y. The Zeitgeber 
in this example is Z(t) — ~(1 + sin(27rf)). 

Remark 1. Here e 'small enough' means that U £ is invertible for which we need that U' £ > 0. This 
depends on the specific form of the Zeitgeber. For the standard Zeitgeber Z(t) = ^(1 + sin(27rt)), small 
enough means e < - . Thus the circle map / M is a diffeomorphism if both e and r\ are smaller than — . ► 

Remark 2. If r\ is not small, but e is small enough so that U £ is invertible, F" 1 is the lift of a circle 
endomorphism, but F^ is multi-valued. This case is similar to case b with the roles of F^ and F^ 1 
interchanged. However in the dynamical system defined with F^ 1 time is running backwards so it 
describes the past rather than the future. Mathematically this is not a problem, but the biological 
interpretation could be problematic. One could force F^ to be single-valued by choosing the smallest 
solution for t = t n+ \ of equation Then F^ again defines a dynamical system, though a discontinuous 
one. ► 

Remark 3. If both e and r\ are not small F^ and F" 1 are multi-valued. In this case we do not have a 
well-defined dynamical system at all. But a similar construction as in the previous remark can be applied 
to define a possibly discontinuous dynamical system. ► 

4 Analysis of the mathematical model 

Here we restrict ourselves to the case that F^ in proposition [2] is the lift of a circle diffeomorphism / M , 
case[a|) on page [5] In that case we can use the rotation number which tells how much on average an initial 
point is rotated along the circle by / M . It is a powerful tool in determining whether has fixed points 
or periodic points. The existence of such points and their dependence on the parameters fi is the main 
topic of this section. Stable fixed or periodic points are the most relevant for our model and we will see 
how they lose stability at certain bifurcations. For background on circle maps and further references to 
the literature see [TJ M EH H3] • 

4.1 Special cases related to the Arnol'd map 

We begin with an example, namely two special cases where F^ can be related to the lift of the Arnol'd 
or standard circle map Aj,a- Here we make a special choice for the Zeitgeber Z, namely Z(t) = \{Y + 
sin(27ri)). The rotation number of the Arnol'd map, depending on the parameters lu and A is well studied, 
so we have quite some information on fixed and periodic points. The (lift of the) Arnol'd map is defined 

as 

Au, a (t) = t + w + A sin(27rt) , 
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for t G R and parameters A, cj e M. In figure [3] we indicate the regions in the (w, A)-plane where the 
rotation number of A w ^ is constant. These regions are called tongues and the general theory of circle 
diffeomorphisms tells us that for parameter values inside the tongues the map has stable fixed or periodic 
points. On the tongue boundaries we have saddle-node bifurcations of fixed points in the main tongue or 
periodic points in the other tongues, see figure [3] 

Recall that the lift of our circle map f^ is given by F^(t) = U n ~ 1 (U E (t + a) — a + T), where U e {t) = t + eZ(t) 
and Z is the periodic forcing. For a special choice of Z and parameters [a, transforms into the Arnol'd 
map by a change of coordinates. We summarize this in the following lemma. 

Lemma 3 (Conjugation to Arnol'd map). Let the periodic Zeitgeber Z be given byZ(t) = |(l+siii(27ri)), 
then 

1. the map F^ £QaT j is conjugate to the Arnol'd map A T+ i £ ^i e , where the conjugation is a rigid 
translation over a, 

2. the map F,^ a r ^ is equal to the Arnol'd map A_ T+ i 1 - j i^. 

Using F^ 1 in the second part of the lemma may seem unnatural in the present context, but if p is the 
rotation number of F , then 1 — p is the rotation number of F. Therefore the second part yields infor- 
mation about the rotation number of Note that the parameter transformation from parameters 
of F to those of A does not involve a. Thus a does not play a role in the bifurcation analysis of these 
two special cases. The proof of the lemma is straightforward and therefore omitted. 




p/q u> 1 w 1 



Figure 3: Left: schematic picture of tongues and hairs. Main tongue emanating at 1, a p : q tongue at 
p/q and a hair at ojq, on the u-axis for the Arnol'd map A Uy \. Right: schematic phase portraits of the 
same family near and on the boundary of the main tongue. The rotation number as a function of u> for 
a fixed value A = Ao is shown in figure [^J 

To analyze these two special cases it suffices to consider the Arnol'd map A u ,\- Let us summarize the 
properties of the latter. If A = the map reduces to a rigid rotation A^fi = Ruj and therefore the rotation 
number is p(A Ui o) = p{Ru) = This is a degenerate situation. But for A ^ the map is no longer 
degenerate. Let us first fix to = ujq and A = A 7^ such that A UOj \ has a rational rotation number |. 
Then A ul0j \ has ^-periodic points. If these points are hyperbolic then there is an open neighborhood of 
(wO)A.o) in the parameter plane such that for all (w, A) in this neighborhood, A u _\ has rotation number 

E 
1 ' 

Let us now consider the line L\ a — {(uj, A ) | Ao 7^ 0, cu 6 [|, |)} in the parameter plane of the Arnol'd 
circle map. From the arguments above it follows that this line segment contains open intervals on which 
the rotation number of the map equals ~. These intervals have to shrink to points when Ao tends to 
zero, because the rotation number of A u ^ equals uj. When we consider the rotation number on the line 
segment L\ as a function of u, its graph, see figure |4j is a so called devil's staircase, see O [23] f° r a 
definition. A further analysis shows that there are saddle-node bifurcations of g-periodic points on the 
boundary of the intervals in L\ , see |24j . Since the map A^^x depends differentiable on the parameters 
there are differentiable pairs of saddle-node curves in the (w, A)-parameter plane emanating from rational 
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points on the line segment L . This forms the structure of tongues in the parameter plane, see figure [3j 
The tongues where the rotation number is | are called p : q-tongues and the tongue where the rotation 
number is 1 (or 0) is called the main tongue. 



P 

1 




UJ 

1 



Figure 4: Devil's staircase: the graph of the rotation number p as a function of U) for the Arnol'd map at 
a fixed value A = Ao. In the main tongue the rotation number is equal to 1. Also see figure^ 

A rational rotation number ^ for the map A Ut \ is constant on closed intervals on the line segment L\ a . 
But the complement of the union of these closed intervals in L\ is not empty, it contains points where 
the rotation number is irrational. Again when Ao tends to zero there are smooth curves with sufficiently 
irrational rotation number ending in irrational points of Lq. These curves are sometimes called hairs. 
The (w, A)-parameter plane of the Arnol'd circle map consists mainly of tongues and hairs, see figure [3j 
The tongues and hairs fill a relatively large region in the parameter plane, when we take a point (u>, A) 
at random, there is a positive probability that it belongs to a tongue, but there is a positive probability 
as well that it lies on a hair. 



4.2 A standard form for the circle diffeomorphism 

The purpose of this section is to show that there is a standard form for every circle diffeomorphism. The 
Arnol'd map for example already is in this form. Identifying the standard form for the map / M allows us 
to conclude that it has a tongues-and-hairs structure in a certain parameter plane for every 1-periodic 
Zeitgeber. Thus the standard form mainly serves a theoretical purpose, for actual computations it is far 
more advantageous to use the expression in equation ^ for the lift of f^. 

Theorem [4] below shows that every lift C of a differentiable circle map of degree one can be written in 
the following form 

C(t) =t + P(t), 

where P is a 1-periodic function. We assume that P is non-zero and non-constant. The Fourier coefficients 
of P may be considered as parameters. Let < P > denote the average of P. Set uj =< P > and 
P = P— < P >. Then uj is the constant term of the Fourier series of P and Pq has zero average. 
Consequently A = max[ 01 ] \Pq\ is non-zero and we may set Poi = jPo- Finally we may write 

C u ,x(t) = t + uj + XP Q1 (t). (5) 

If we now interpret uj and A as parameters, then there is a tongues-and-hairs structure in the (uj, A)-planc. 
Note that we recover the Arnol'd map by setting Poi(i) = sin(27rf). In the more general family of equation 
([5]) the tongues may have a richer structure than those of the Arnol'd map. For example there may be 
more saddle-node curves inside a tongue, as we will see in the next section. For a detailed description of 
such phenomena see [5]. 

Recall that F^(t) = U~ 1 (U E (t + a) — a + r), where U e (t) =t + eZ(t) is a lift of the circle diffeomorphism 
f^. First we give a standard form for the map F^. 
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Theorem 4 (Standard form). Suppose that r\ is small enough so that U~ x exists. Let v = (a, ft, a,r) be 
new coordinates in parameter space with e = crcos/3 and rj = a sin (3. Then there are smooth functions u> 
and X of the parameters v with u(y)\ a= Q = r, \{v) = a and a 1-periodic smooth function i?„ with zero 
average, smoothly depending on parameters v such that 

F^(t) = t + uj(v) + \(v)R„(t). 

The theorem shows that after a transformation of parameters, takes the standard form of a circle 
map. In particular we have the result that there is a tongues-and-hairs structure in the (r, er)-plane, for 
each value of (3 e [0, f ) and each value of a £ [0, t). 

The function R depends on the Zeitgeber Z, but also on all parameters v. The latter will appear in the 
coefficients of the Fourier series of R. 

Remark 4. Setting rj = in fi = (e,r/,a,T) corresponds to setting f3 = in v = (a, f3,a,r). Let 
<Z>= f Q Zit) dt be the average of Z then wc have 

F( E ,o,a,r) =t + r + eZ(t + a) =t + r + e <Z> +e(Z(t + a)- <Z>). 

If we also take Z(t) = |(1 + sin(2-7rf) we recover the result of lemma [SJ ► 

Remark 5. We expect that for most values of j3 the tongues in the (r, cr)-parameter plane are non- 
degenerate, that is they intersect transversely at the vertices. We even expect this for Z(t) = |(1 + 
sin(27ri)), since in the standard form of F, R v will have a Fourier series rather than a Fourier polynomial. 
Also see [H]. ► 

Let us now look at the position of the main tongue. We assume that the Zeitgeber has the following 
Fourier series 

Z(t) =c + ci sin(27rf) + °k sin(2vr(H + j k )), (6) 

k>l 

with coefficients Ck G K and 7^ s [0, 1]. Note that by a time shift we can always achieve that 71 = 0. 
Then we have the next result. 

Proposition 5 (Boundaries of main tongue). Let the Zeitgeber be as in equation^ then the boundaries 
of the main tongue of the map F^^ of are given by 



t± = 1 — a 



c (cos - sin (3) =p c 1 y/l - cos(2a7r) sin(2/3)j + 0(a 2 ). (7) 



For our standard example of a Zeitgeber Z(t) = \(1 + sin(27rf)) the boundaries of the main tongue are 
given by equation ^ with cq = c\ = | and without the 0(a 2 ) term because all other coefficients are 
equal to zero. 



4.3 Fixed points of the diffeomorphism f ^ 

Existence of stable fixed points of the circle diffeomorphism f ^ is one of the main questions. Recall that 
such points correspond to entrainment in the biological model. Therefore we take a closer look at such 
points. From the previous sections we know that fixed points exist for parameter values in the main 
tongue. They are easily characterized as follows, t = t* is a fixed point of /„ if / M (i*) = t*. However, as 
noted before, for practical computations it is more convenient to use a lift of / M . Using the notation 
from proposition [2] a point t — t* is a fixed point of if 

F^{t*) = t* + 1, 

which is equivalent to 



U e (t* +a)-a + T = U v (t* + 1). 



(8) 
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Even if e and rj are small enough so that U e and U„ are invertible, this equation has several solution 
branches. Therefore we can not solve ([8| explicitly for t* . Since we are mostly interested in the dependence 
of t* on r we shall content ourselves with 

t = U n (t*) - U £ (t* + a ) + a + l. (9) 

Thus for fixed values of e, rj and a we obtain t as a function of t* . 

Now that we have characterized the fixed points of let us determine their stability. The fixed point 
t = t* is stable if \F'(t*)\ < 1. After a short computation we find 

F>(n = u ^+ a \ do) 

From equations ^ and ( 10 ) we almost immediately obtain an alternative characterization of stability of 
the fixed point t* . Two examples of application of the following lemma are given in figure [5] 

Lemma 6 (Stability). On solution branches of equation (J8|, the fixed point t* is stable or unstable when 
t* as a function of r is increasing or decreasing. 

Let us again consider our standard Zeitgeber Z(t) = |(1 + sin(27rf)). For fixed e, 77 and a there are two 
values t* and t* 2 for which F'^(t*)\ = 1 and |-F^(t*)| < 1 if t* € [t*,^]. We recover the boundaries of the 
main tongue by setting n = U v (t*) — U E (t* + a) + a + 1, then a stable and an unstable fixed point of 
the map / M exist for r on the interval [ti,T2]. Motivated by this example one could conjecture that for 
every 1-periodic Zeitgeber with only two extrema in one period there are precisely two fixed points for 
parameter values in the main tongue. This turns out to be true when an extra condition is imposed. 

Proposition 7 (Number of fixed points). Let Z be a 1-periodic Zeitgeber with one maximum and one 
minimum on the interval (0, 1) . Then the circle map f ^ has precisely two fixed points for parameter values 
in the main tongue if (Z") 2 — Z' ■ Z'" does not change sign. 

For our standard example of a Zeitgeber Z(t) = |(l + sin(27rf)) the quantity (Z") 2 — Z' ■ Z'" is equal to 1. 
The Zeitgeber Z(t) = |(| + sin(27rf) + ± sin(47rt)) also has two extrema, but the quantity {Z") 2 - Z' ■ Z'" 
changes sign on [0, 1]. See the appendix for further implications. With a general periodic Zeitgeber there 
may be more than one stable fixed point for parameter values in the main tongue. This occurs in general 
when the Fourier series of the Zeitgeber contains more than just one term like in our standard example. 
Let us look at a Zeitgeber with the following Fourier polynomial Z(t) = ^(2 + sin(27rf) + cos(47rf)), then 
there are four solutions of F'(t*)\ = 1. Let us sort them such that r, = U v (t*) — U £ {t* + a) + a + 1 
increases with i. Then for fixed values of e, rj and a there are at most four fixed points of for each 
re [n, T4] . We graphically represent the results in figure [5j 



4.4 Examples of tongues-and-hairs for different Zeitgebers 

Here we collect some examples of tongues-and-hairs figures showing differences and similarities with the 
prototype figure of tongues in the Arnol'd circle map. Therefore we start with the latter, see figure [6j In 
this family the tongues are fixed, but in our family f^ , the tongues in the (t, <r)-plane still depend on the 
values of j3 and a. Moreover they also depend on the Zeitgeber. In order to keep the number of pictures 
limited we only show tongues for two different Zeitgebers, see figures [7] and [8] As is to be expected form 
the existence of a standard form, see section |4.2| the pictures are qualitatively the same. That is near 
the r-axis or w-axis in the Arnol'd map. The width and the growth of the width of the tongues depends 
on parameters in the map. This is shown in figure [9j in the next section, where its biological relevance is 
discussed. 
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Fi gure 5: Fixed points t* of the map as a function of t . The solid cuvve vepvesents a stable fixed point, 
the dashed curve an unstable one. At parameter values t = Ti there are saddle-node bifurcations, also see 
figure^^ Parameters e, n and a are fixed. Left: Zeitgeber is Z(t) = ^(1 + sin(27rf)). Right: Zeitgeber is 
Z(t) = ^(2 + sin(2vrf) + cos(47ri)). 







1.0 



u 1.5 



Figure 6: Tongues for the Arnol'd map A Uj \ in the (a;, \)-plane. The algorithm to compute the pictures is 
based on the rotation number. The latter is defined for sufficiently small values of A only. Therefore the 
tongue boundaries in the picture are not well defined for large values of X. The effect of this phenomenon 
is even more visible in the following picture. 




0.5 



1.0 



1.5 



Figure 7: Tongues for the circle map f^ in the (r, a) -plane, with standard Zeitgeber Z(t) — i(l+sin(27rf)). 
The values of the parameters /1 = (e, 77, a, r) are e — a cos (3, rj — a sin /3 with j3 — | and a — 0.3. 
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5 Discussion and conclusion 



Model 

The biological model for a single pacer cell with a periodic Zeitgeber as described in section [2] leads to 
quantitative mathematical model, namely a dynamical system, which we now summarize. The heart of 
the mathematical model is a circle map / M that yields once to is given, a sequence of relative times t n 
corresponding to the beginning of an activity interval of the pacer cell. The sequence to, t\ = f^to), t% = 
fij,(t\), . . . is called the itinerary of to- In this setting time is relative to the period T of the periodic 
Zeitgeber that acts as a stimulus to the pacer cell. It is convenient to use T as a unit of time, that is 
we scale time so that T = 1. The circle map / M depends on parameters fi — (e,r), a,r), where r is the 
intrinsic period of the pacer cell, measured in time unit T = 1. Furthermore e and r\ determine phase 
delay and phase advance and a determines the length of the intrinsic activity interval, see section [2] for 
an explanation of these terms. 



Dynamics 

In the description of the dynamical behavior of the map / M there is a difference between / M being 
invertible (diffeomorphic) or not. We mainly restrict to the invertible case. Furthermore we restrict to 
typical dynamics, see [9] for a precise definition. In the present case this comprises fixed points, periodic 
points and quasi-periodic points. Note that the itinerary of the latter consists of an infinite sequence. 

The simplest kind of dynamics is a fixed point of which means that the activity interval of the pacer cell 
always starts at the same relative time. The existence of such points implies the possibility of entrainment. 
The next simplest kind of dynamics is a periodic point of f^. This means that there is a to such that in 
the sequence to, t± = / M (£o), ■ ■ ■ ,t q — / M (i<j-i) the last point t q is again equal to to for some fixed q. In 
general these q onsets of activity occur in p periods of the Zeitgeber. We call this synchronization of the 
pacer cell, a generalization of entrainment. The periodic point to is called p : q periodic. The last kind 
of typical dynamics for the map , is quasi periodicity. The point to is quasi-periodic if the itinerary of 
to densely fills the circle or interval [0, 1] depending on how we represent the map. 

The kind of dynamics the map / M exhibits depends on the values of the parameters. All this is nicely 
organized in parameter space in wedge shaped regions called tongues, one for each pair (p,q), see figure 
[3j However, here we need the restriction that / M is invertible. For parameter values in the (p, q) tongue, 
the typical dynamics of is p : g-periodicity. A special role is played by the tongue for (p, q) = (1,1) 
called the main tongue. Fixed points are the typical dynamics of / M for parameter values in the main 



tongue. At the boundaries of the tongues the map / M has a saddle- node bifurcation, see section 4.1 For 



parameter values outside the tongues the dynamics of the map is quasi-periodic, this occurs on hairs 
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in the (r, cr)-parameter plane. If /„ is not invertible, tongues still exist but generally overlap so that 
coexistence of periodic points with different periods becomes possible. 



Entrainment 



For parameter values in the main tongue, the tongue with (p, q) — (1, 1), the map /„ has fixed points. 
There are at least two such points of which one is stable and the other is unstable. The stable one 
corresponds to entrainment of the pacer cell. Let us discuss the shape of the main tongue in some more 
detail to find the parameter values for which entrainment occurs. 

The parameter space is four dimensional and in it tongue boundaries are hyper-surfaces. The situation 
becomes simpler if we do not use coordinates (e, 77, a, t) but (<r, /3, a, r) with e = a cos {3 and rj = a sin j3. 
Then a satisfies a 2 = e 2 + rj 2 and measures how strongly the Zeitgeber stimulates the pacer cell, while 
/? determines the ratio of e and 77. In these coordinates in parameter space the boundaries of the main 
tongue are given by equation Q. As we can see from the standard form in theorem|4]the main parameters 
are r and a. In the (r, <r)-plane we find the tongues, see figure [9j whose detailed shape depends on a and 
/?. The range of entrainment is given by the interval (r_, r + ), with t± given in equation ([7]), when other 
parameters are kept fixed. There are many biological experiments/observations supporting the existence 
of bounded ranges of entrainment, see |10| . In practice one cannot vary the intrinsic period r of the 
pacer cell. But varying the period of the Zeitgeber T at a fixed value of r is equivalent to varying r 
and fixing T in the model, see 5, 30] [34]. However, biological evidence exists that the intrinsic period 
varies among individual pacer cells while they can still entrain to a Zeitgeber with a 24-hour period, see 

[H mini [31 H2 • 



7T 




T_ 1 T + T T_ 1 T + T 



Figure 9: Left: the main tongue of the circle map f u> where v = (tr, j3, a, r) with a = 0.3 and f3 — /3q fixed. 
For o = o~o the range of entrainment is the width of the tongue at a = o~q, namely the interval (r_,T + ). 
Right: the range of entrainment (t_,t + ) for fixed a — a and f3 varying in (0, ^). The parameter (3 
determines the ratio of e and rj since e = trcos/3 and rj = crsin/3. 



In section 4.3 on fixed points of the map we noted that the position of stable fixed points is an 



increasing function of the intrinsic period r. This has been observed by various authors, for example 



[2|I3T] and for particular organisms by [U O [HI [30] . In figure 10 we show data from [20] essentially giving 
the relative onset times of the activity interval as a function of the intrinsic period, of several mutants of 
the fungus Neurospora crassa. 

Both position and length of the interval (t_,t+) depend on a and /?. As we see from figure [9] the range 
of entrainment is in general not centered at r = 1. Here f3 is the most important parameter. If /? < ? or 
equivalently e > rj, then phase delay is larger than phase advance and the range of entrainment is shifted 
towards intrinsic periods r smaller than the period T of the Zeitgeber. If /3 > ^, the range of entrainment 
is shifted in the direction of intrinsic periods r larger than the period T of the Zeitgeber. The model 
has two extreme cases, one for /3 = or equivalently rj = (only phase delay), where only pacer cells 
with intrinsic period r less than the period T of the Zeitgeber can be entrained. The other one is for 
P = j or equivalently e = (only phase advance), where r must be larger than T for entrainment to 
occur. A phenomenon related to this skewness has been observed in several nocturnal rodents |11] where 
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Figure 10: Relative times of onset of activity interval for four mutants of Neurospora crassa. The genetic 
types are indicated by markers, see J3U\j for an explanation of frql, frq+, frq7 and frq9&10. 

a decreasing center of the range of entrainment corresponds to increasing phase delay and decreasing 
phase advance. 

As mentioned before, in the main tongue, the map / M has at least a pair of fixed points, one stable and 
one unstable. However more pairs may exist, leading to curves of saddle-node bifurcations inside the 
main tongue, see figure [5j Thus upon varying r a fixed point may lose its stability and the system jumps 
to another stable fixed point. As long as parameter values remain in the main tongue this is the only 
possibility. Another possibility is to keep the parameters fj, fixed but vary a parameter in the Zeitgeber. 
For a biological example possibly related to such a mechanism see [35]. 

Synchronization 

Apart from entrainment, the model also shows the possibility of the more general phenomenon called 
synchronization. This occurs for parameter values in the p : q tongues in the (r, er)-planc. Then the circle 
map ffj, has a g-periodic orbit consisting of points to,ti,... , t g _i indicating the beginnings of q activity 
intervals of the pacer cell in p periods of the Zeitgeber. In the p : q tongues we have the same phenomena 
as in the main tongue, the only difference is that they apply to periodic points instead of fixed points. 

There is an example of 2:1 periodic point, the fungus Neurospora sp has one activity interval in two 
periods of the Zeitgeber, see |25j . Here the Zeitgeber is a temperature stimulus. 

No synchronization 

For parameter values outside the tongues the circle map has quasi-periodic orbits. This corresponds to 
quasi-periodic occurrence of activity intervals of the pacer cell. In practice such behavior may be hard to 
distinguish from periodic behavior with a long period. 

Future directions 

There are several ways to generalize or extend the current model for a single pacer cell with a periodic 
Zeitgeber. We first concentrate on the biological model of section [2] in view of our future goal to describe 
a collection of interacting pacer cells, stimulated by a Zeitgeber. 
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Single pacer cell, circle map In the present analysis we restricted to the case that the map / M is a 
circle diffeomorphism (invertible map). Then the (t, er)-plane is divided into tongues with well-defined 
periodic dynamics and hairs with quasi-periodic dynamics. This only occurs when the stimulus of the 
Zeitgeber on the pacer cell is relatively weak. Allowing a stronger stimulus, the map becomes an 
endomorphism (non-invertible map) with far richer bifurcation scenarios, see for example [8]. However, 
it is not clear whether the subtleties of the endomorphism case are in accordance with the coarseness of 
the underlying biological model. Therefore we restrict to the simplest properties of each model. 

Single pacer cell, torus map Thus it seems more fruitful to generalize in another direction and 
consider a quasi-periodic Zeitgeber. This will lead to a torus map instead of a circle map. Such a 
generalization is also more relevant for our aim to study a collection of pacer cells, by first considering a 
single cell with an asymmetric interaction with its environment. The latter stimulates the pacer cell, but 
not vice versa. 

Collection of pacer cells, no interaction In a model for a collection of pacer cells we can already 
use the results for a single cell. As a first model let us assume that the cells are stimulated by an 
external periodic Zeitgeber but do not interact. Although biologically not particularly relevant, it is a 
step in gradually sophisticating the model. Furthermore suppose that there are n cells, characterized by 
parameter values (£j, rji, oti, Ti) in the main tongue for i — 1, . . . , n. That is we apply the current single 
pacer cell model for each cell. Then the collection will be entrained to the Zeitgeber albeit that each cell 
has its own onset time and length of activity interval. The collective behavior though, will be periodic. 

Collection of pacer cells, with interaction We conjecture that, starting with the previous model, 
for a sufficiently weak interaction there will still be entrainment. Nevertheless it will be interesting to 
consider a model for interacting pacer cells without Zeitgeber as well. Here we may take the pacer cells 
nearly identical in the sense of the previous paragraph. However there is biological evidence that there 
are different types of pacer cells [THl |3B] . In further extensions one could again include a periodic 
Zeitgeber. Since in a model for a collection of pacer cells the Zeitgeber acts solely as an external stimulus, 
it seems most natural to restrict to periodic Zeitgebers. However we may wish to include both daily and 
seasonally variations. This would imply that the period of the Zeitgeber is a year. Another possibility is 
to stick to Zeitgebers with a 24-hour period and use methods for slowly varying parameters to include 
seasonal changes. 

A Bifurcations 

The tongues in the (<t, r)-parameter plane are determined by saddle-node bifurcations. But there may also 
be other bifurcations even when the map f ^ is a diffeomorphism, in other words invertible. The reason 
is that we have many parameters. Considering a as a relatively unaccessible parameter and keeping it 
fixed we may still vary er, t and j3. Then we have 3-dimensional tongues in (cr, t, /3) -parameter space. It 
is an almost straightforward consequence of theorem [7] that there are curves of pitchfork bifurcations in 
this three dimensional parameter space, emanating from rational points on the r-axis. 

Corollary 8 (Pitchfork bifurcation). Let Z be a 1-periodic Zeitgeber. If (Z") 2 — Z' ■ Z'" has a simple 
zero, then parameter values exist for which / M has a pitchfork bifurcation. 

Proof. From the proof of proposition [7] we see that the number of solution branches of equation Q does 
change if Z"(t)' 2 - Z'(t) ■ Z"'(t) has a simple zero. ■ 

However, (3 may be considered as an unaccessible parameter as well. But the Zeitgeber may also depend 
on a parameter. We may in particular view seasonal change, which is slow compared to the 24-hour 
period, as parameter dependence. The quantity (Z") 2 — Z' ■ Z'" may change sign depending on this 
parameter, so we find again pitchfork bifurcations. 
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If the map f ^ is not a diffeomorphism (not invertible) then there are numerous other bifurcations, see 
[S]. This happens for relatively large values of a. On varying r for fixed values of the other parameters 
in the main tongue, one generally finds a number of period doublings followed by the same number of 
period halvings (or in opposite order). The reason is that there are curves of period doublings in the 
(r, (r)-parameter plane with local minima, considered as functions of r, that are transversally crossed. 
For more details we refer again to [8]. An example of this phenomenon is shown in figure 

0.5: 



t 



11 





0.72 

Figure 11: Period doublings and halvings for the map f ^ with Zeitgeber Z(t) 
Shown are the positions of periodic points as a function of the intrinsic 
a = 0.5, e = 0.4 and rj = 0.17. 



t 1.05 

•sin(27rf) + 2cos(47rf). 
Parameter values are 



B Proofs 



Proof of proposition [1| The main point we have to show is that equation |2]) can be solved uniquely with 
respect to t. Using the expression for in (JlJ we obtain after some rearranging 



t n+1 + rjZ(t n+l ) =t n + a + eZ(t n + a) - a + r. 



(11) 



Since Z is 1-periodic, U e has the property U e (t + 1) = U E (t) + 1, for all t. Then the equation for t n+1 
reads 

U v (t n +i) = U e (t n + a)-a + r. 

Introducing the operator T a which takes a function / into T a f — T_ a o f oT a where T a is just translation 
over a, that is T a (t) = t + a we can write the equation for t n+ i as 



U v (t n+1 ) = T a U E (t n ) + t. 



(12) 



Solvability of (12) depends on the value of rj. If f] is small enough, is invertible and we write 

t n+ i = F^{t n ) = U-\T a U £ (t n ) + t). 

Then F M is a differentiable map with the property F^(t + 1) = F^(t) + 1, which means that is the lift 
of a circle map of degree one. Note that F^ depends on the parameters fi = (e, 77, a, t). I 

Proof of theorem^ A consequence of F^it + 1) = F^(t) + 1 is that F^(t) — t is 1-periodic, thus there is 
a 1-periodic C°° function P v such that F^ (t) = t + P v (t) . P v has a Fourier series so we may split off the 

constant term and we write P v (t) = u{v) + Q u (t) with u)(v) = f Q P v (t) dt, then uj is a C°° function of v. 
Furthermore Q v (t) — P v (t) — u)(v) so that Q v is a 1-periodic C°° function with Q v {t) dt = 0. So far 
F/j,(t) = t + uo{v) + Q v {t). For a = we have F M (t) = t + r, so u(Q, (3, a,r) = r and Qro,^,a,r)(*) = 0. Using 
the division property of C°° functions, a 1-periodic, C°° function R v exists such that Q v {t) = <jR u {t). 
Finally F {t7COS ^ sin 0, a , T ){t) = t + u(v) + crR v {t). ■ 
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Proof of proposition^ Recall that fi = (s,r),a,r) and v — (er, (3, a, r) with e = a cos /3 and 77 = crsin/3. 
Then for small a and assuming that 1 — r = 0(c) we get 

F„ {v) = U_ v {U e (t + a)-a + T)+0(a 2 ) 

= t + r + eZ{t + a) - j]Z{t + t) + C>0 2 ) 
= t + r + eZ{t + a) - t]Z(t) + 0(cr 2 ). 

The Zeitgeber Z has the following Fourier series Z(t) = c Q + ci sin(27rf) + Ylk>i ° k sin(27r(fci + 7^). After 
a near identity transformation followed by a time shift we obtain that F^^ is equivalent to 

G^{u){t) =t + T+(e~ n)c + ec x sin(27r(t + a)) - ija sin(27ri) + C(cr 2 ) 

It almost immediately follows that the boundaries of the main tongue of G and thus of F are as stated 
in the lemma. I 

Proof of proposition [7| The number of fixed points may change if the number of solution branches of 
equation (J8J): U e {t + a) — a + r = U v (t + 1) changes. This is equivalent to a changing number of extrema 
of r = U n (t) - U s (t + a) + a + 1. Using U e (t) = t + eZ{t) we get r = i]Z(t) - eZ{t + a) + a + 1. Thus 
the number of extrema of r changes at parameter values for which 

r,Z'{t)-eZ'{t + a) = 
r)Z"(t) -eZ"(t + a) = ' 

This equation has trivial solutions for e and 77 only if Z'(t) ■ Z"(t + a) — Z'(t + a) ■ Z"(t) ^ 0. We rewrite 
this as 

Z'{t) Z'(t + a) 



Z"{t) ' Z"{t + a) 

,AA — 

Z"(t) 



Z' (t) 

This inequality holds if h(t) — z >>(t\ ^ s injective. Therefore a sufficient condition is that h does not 



change sign. From 

_ Z"(tf-Z'(t)-Z"'(t) 

[, ~ z»{ty 

the result follows. 
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